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Abstract — A key issue in physical-layer network coding (PNC) 
is how to deal with the asynchrony between signals transmitted by 
multiple transmitters. That is, symbols transmitted by different 
transmitters could arrive at the receiver with symbol misalign- 
ment as well as relative carrier-phase offset. In this paper, 1) 
we propose and investigate a general framework based on belief 
propagation (BP) that can effectively deal with symbol and phase 
asynchronies; 2) we show that for BPSK and QPSK modulations, 
our BP method can significantly reduce the SNR penalty due to 
asynchrony compared with prior methods; 3) we find that symbol 
misalignment makes the system performance less sensitive and 
more robust against carrier-phase offset. Observation 3) has the 
following practical implication. It is relatively easier to control 
symbol timing than carrier-phase offset. Our results indicate that 
if we could control the symbol offset in PNC, it would actually 
be advantageous to deliberately introduce symbol misalignment 
to desensitize the system to phase offset. 

Index Terms — physical-layer network coding, network coding, 
symbol synchronization, phase synchronization 

I. Introduction 

Physical-layer network coding (PNC), first proposed in |IT], 
is a subfield of network coding I2j that is attracting much 
attention recently. The simplest system in which PNC can be 
applied is the two-way relay channel (TWRC), where two end 
nodes exchange information with the help of a relay node in 
the middle, as illustrated in Fig. [TJ^)- 

Compared with the conventional relay system, PNC doubles 
the throughput of TWRC by reducing the needed time slots 
from four to two. In PNC, in the first time slot, the two end 
nodes send signals simultaneously to the relay; in the second 
phase, the relay processes the superimposed signals and maps 
them to network-coded information for broadcast back to the 
end nodes. From the network-coded information, each end node 
then makes use of its self information to extract the information 
from the other end node. 

A key issue in PNC is how to deal with the asynchrony 
between signals transmitted by multiple transmitters. That is, 
symbols transmitted by different transmitters could arrive at the 
receiver with symbol misalignment as well as relative carrier- 
phase offset. 

Many previous works {e.g., Ill, |[3l, ID) found that symbol 
misalignment and carrier-phase offset will result in appreciable 
performance penalties. For BPSK modulation, IT] showed 
that the BER performance penalties due to the carrier phase 
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asynchrony and symbol offset are both 3 dB in the worst case. 
For QPSK modulation, the penalty can be as large as 6 dB in 
the worst case when the carrier-phase offset is tt/A @|. 

These earlier investigations led to a common understanding 
that near-perfect symbol and carrier-phase synchronization are 
important for good performance in PNC. We believe that this 
is a misconception, and that asynchronous PNC can have good 
performance if we use the right methods to deal with the 
asynchrony. 

The study of BPSK PNC in [B, 0, for example, made 
use of suboptimal decoding methods at the relay for the asyn- 
chronous PNC case. Furthermore, the joint effect of symbol 
and phase asynchronies was not investigated. In this paper, 
we propose an optimal maximum-likelihood (ML) decoding 
method making use of a belief propagation (BP) algorithm that 
addresses symbol and phase asynchronies jointly. We find that 
the BER performance penalty can be reduced from 3 dB (using 
the method in HI, 13) to less than 0.5 dB (using our method 
here) for the worst case. 

In QPSK PNC, the penalty is larger than 6 dB |f4| only when 
the symbols are perfectly aligned and when the phase offset is 
7r/4. We find that using our optimal BP decoding algorithm, 
when there is a 1/2 symbol misalignment, the penalty is less 
than 1 dB (compared with the case when symbol and phase are 
perfectly synchronized). Additionally, an interesting result is 
that with the 1/2 symbol misalignment, the spread of penalties 
under various carrier-phase offsets is no more than 0.5 dB. This 
means that symbol misalignment has the effect of desensitizing 
the performance of the system to carrier-phase offset. The 
practical implication is as follows: 

In practice, it is easier to control symbol timing than carrier- 
phase offset because of the relative time scales involved. If we 
could control the symbol offset in PNC {e.g., (S) presented 
a method to synchronize symbols from different sources), 
it would actually be advantageous to deliberately introduce 
symbol misalignment so that the system is robust against 
uncontrollable phase offsets. 

The remainder of this paper is organized as follows: Section 

II introduces the system model adopted by this paper Section 

III classifies and defines synchronous PNC and asynchronous 
PNC. Section IV presents our BP ML decoding method for 
dealing with asynchrony in PNC. Numerical results are given 
in Section V. Finally, Section VI concludes this paper 

II. System Model 

We study the two-way relay channel as shown in Fig. [T] in 
which nodes A and B exchange information with the help of 
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Fig. 1. System model for two way relay channel. 

relay node R. We assume that all nodes are half-duplex, i.e., 
a node cannot receive and transmit simultaneously. We also 
assume that there is no direct link between nodes A and B. 

We consider a two-phase transmission scheme consisting of 
an uplink phase and a downlink phase. In the uplink phase, 
nodes A and B transmit packets to node R simultaneously. In the 
downlink phase, based on the overlapped signals received from 
A and B, R constructs a network-coded packet and broadcast the 
packet to A and B. Upon receiving the network-coded packet, 
A (B) then attempts to recover the original packet transmitted 
by B (A) in the uplink phase using self-information 11 1. 

Consider the uplink phase again, if only A transmitted, then 
the received complex baseband signal at R (i.e., the received 
signal after down-conversion from the carrier frequency and 
low-pass filtering) would be 



N 



VB{t) = 22 hAXA[n\p{t - nT) + WB.{t), 



where Ha = VPa is the received signal amplitude; 
{xA[n])n=i,....N are the symbols in the packet of A; p{t—nT) is 
the pulse shaping function for the baseband signal; and wj^ (t) 
is Gaussian noise with double-sided power spectral density 
Sw{f) = No/2 for both the real and imaginary components 
within the baseband of interest. In the PNC set-up, A and 
B transmit simultaneously. In this case, the received complex 
baseband signal at R is 



N 



VR 



[t) = ^ {hAXA[n]p{t — nT) + hBXB[n]p{t — A — nT)} 



WR{t), 



(2) 



where Hb = y/PBe^"^ {(f> is the relative phase offset between 
the signals from A and B due to phase asynchrony in their 
carrier-frequency oscillators, and difference in the delays of 
the two uplink channels); {xB[n])n=i n are the symbols in 
the packet of B; and A is a time offset between the arrivals 
of the signals from A and B. Without loss of generality, we 
assume the signal of A arrives earlier than B. Furthermore, we 
assume A is within one symbol period T. Thus, < A < 
rOl We refer to A and <f) as the symbol and phase offsets (or 

'if A is more than one symbol period, we could generalize our treatment 
here so that A'^ is larger than the number of symbols in a packet. The 
packets will only be partially overlapping, with non-overlapping symbols at 
the front end and tail end. Essentially, our assumption of A being within 
one symbol period implies that we are looking at the "worst case" with 
maximum overlapping between the two packets. When there are additional 
non-overlapping symbols at the front and tail ends, the decoding will have 
better error probability performance. 



misalignments) at R, respectively. Note that we have assumed a 
flat fading channel, and the channel gains stay constant within 
a packet duration. For simplicity, we assume power control so 
that Pa = Pb = P- Furthermore, for convenience, we assume 
time is expressed in unit of symbol duration, so that T = 1. 
Then, we can rewrite (|2) as 

N 

ynit) =^Y1 {^A[n]p{t ~ n) + XB[n]p{t - A - n)e^^} 

Tl=l 

+ WRit). (3) 

In general, the pulse shaping function p{t) can take different 
forms. To bring out the essence of our results in the simplest 
manner, throughout this paper, we assume the rectangular pulse 
shape: p{t) = rect{t) = u{t + 1) — u{t). 

A critical design issue is how node R makes use of 
to construct a packet for broadcast to nodes A and B in 
the downlink phase. In this paper, we assume that R first 
oversamples yR{t)/^/P to obtain 2N + 1 signal samples. The 
2N + 1 signal samples will be used to construct an A^-symbol 
packet for broadcast to A and B. The oversampling procedure 
is similar to that in 16] and described below. For n — 1, N, 



I r{n-l)+A 

yR[2n - 1] = -r / {xA[n\ 

^ J (71-1) 



{xA[n]+XB[n- - 
ly'l' + WR[2n~l], 
{xA[n]+ XbMc^''' 4 
XA[n\ + XBinle^"^ + WR[2n\, 



WRit) 



l{n- 

XA[n] + XB[n 
1 



)dt 



1- A 



WR{t) 



(ri-l)+A 



P 

(4) 



)dt 



(1) and 



yR[2N + i] 



N+A 



(^XBiNy'f' + WR{t)/Vp)dt 

1], 



AVPjn 
= XB[Ny>' + wr[2N ■ 

where xb[0] ~ 0, and WR[2n — 1] (also wr[2N + 1]) and 
i(;fl [2n] are a zero-mean complex Gaussian noise with variance 
No/i2PA) and No/ (2F(1 - A)), respectively, for both the 
real and imaginary components. 

Fig- nib) is a schematic diagram of our system. This paper 
adopts the following notation: 

» Xi = {xi[l], Xi[2], a;i[Af]) denotes the source packet of 

node /, i G {A, B}; 
. YR = {yR[l],yR[2],...,yR[N],yR[N+l],...,yR[2N + l]) 
denotes the received packet (with the aforementioned 
oversampling) at relay node R; 
. Wr ^iwR,[l],WR[2],..., WR [N] ,WR,[N + 1],..., WR [2N + 1] ) 
denotes the receiver noise at node R; 

• Xr = (a;7{[l], a:;i^[2], Xfl[Af]) denotes the network- 
coded packet at relay node R; 

' Yi = {yim,yi[% -^yAN]) denotes the received PNC 

packet at node i G {A, B}\ 
» Wi = {wi[l],Wi[2], ...,Wi[N]) denotes the receiver noise 

at node /, i £ {A, B}; 

• Xi ~ {x.i[l], x.i[2], ...,Xi[N]) denotes the decoded source 
packet of node ;, i G {A, B}, at the other end node; 

where N is the number of source symbols. 
In the above, for BPSK, a;i[-], a;i{[-] G 
{-1,1}; for QPSK, Xi[-],Xi[-],XR[-] G 
{(1 + j)/ A (-1 + j)/V2,{-l - j)/V2Xl - j)/V2} . 
In addition, yA-]-,yR[-],wA-],WR[-] G C. 
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III. Synchronous vs. Asynchronous PNC 

A. Classification 

Table H] shows the four possible cases in PNC systems. Case 
1 is the perfectly synchronized case studied in JT]; Case 2, in 
0; and Case 3, in fS), Q, lU. This paper proposes a general 
scheme to tackle all four cases under one framework. 

TABLE I 
Four cases of PNC systems 





A = 


A ^0 


(/) = 


(/, = 0, A = (Case 1) 


= 0, A 7^ (Case 2) 


ct>^0 


^ 0, A = (Case 3) 


^ 0, A 7^ (Case 4) 



B. Review of Synchronous PNC ( Case 1 in I///-AI ) 

Definition 1 (Synchronous PNC): Synchronous PNC is the 
case where the two end nodes transmit their packets Xa and 
Xb in a synchronous manner so that the relay node R receives 
the combined signals with = and A = 0. The received 
baseband packet at R'lsYn — Xa+Xb + Wr with symbols. 
Node R then transforms Yji into a network-coded packet X^ — 
KXr) for transmission in the downlink phase. 

For this case, with reference to (@}, since A = 0, the variance 
of the noise term Wfl;[2n — 1] is infinite, and the signal is 
contained only in the even terms yft\2n]. Thus, we can write 

yR[2n] = XA[n] + XB[n] + wii[2n], (5) 

where n = \, ...,N , and i(j7^[2n] is zero-mean Gaussian noise 
with variance = for both the real and imaginary 

components. 

For BPSK, Xi[n] € {—1,1}. Only the real component of 
Wfl[27i] needs to be considered. For QPSK, since we are con- 
sidering a synchronous system, the in-phase and quadrature- 
phase components in (|5]l are independent; it can therefore be 
considered as two parallel BPSK systems. In the following, we 
assume BPSK. 

Let us consider a particular «, and omit the index n in 
our notation for simplicity. The a posteriori probability of the 
combination {xa, xb) is given by 

'Py{VR\xA,XB) 



Vt:{xa,xb\vr) 
1 



4Pr(2/i?) 



: exp 



XA - Xb)] 



(6) 



4Pr(yfl)V27rcr2 "'^ 2(t2 

Let us use Xi — Ito represent bit and Xi = —1 to represent 
bit 1. Suppose that the downlink transmission also uses BPSK. 
For PNC output, xr, let us consider the XOR mapping: xr = 
XA © Xb ID- Then, xr = 1 if = xb, and xr = —1 if 
xa 7^ Xb- The following decision rule is used to map yR to 
Xr: 

xr = 1 

1, XB = l|y_R) + P^ixA = -1, a;_B = -Mvr) 



Pr{xA 
PiixA 
exp 



1,XB = -IIvb) + Pt{xa = -l,XB = IIvr) 



{yn ~ 2)^ 
2^2 



exp 



2cr2 



Xr = -1 



exp 



Xr = ~1 

yjL 

2^2 
(7) 



C. Asynchronous PNC ( Case 2, 3 and 4 in I///-AI I 

Definition 2 (Asynchronous PNC): Asynchronous PNC is 
the case where the two end nodes transmit their packets Xa 



and Xb in an asynchronous manner so that the relay node 
R receives the combined signals with 7^ or A 7^ 0. The 
received baseband packet at R is Yr = Xa + Xr + Wr with 
2N + 1 symbols. Node R then transforms Yr into a network- 
coded packet Xr ~ /{Yr) for transmission in the downlink 
phase. 

The asynchronous case has been studied previously. How- 
ever, suboptimal PNC mappings Xr = KXr) were consid- 
ered. Refs. IT] and [S) argued that the largest BER performance 
penalty is 3 dB for BPSK modulation (for both phase and 
symbol asynchronies). These results, however, are based on 
suboptimal decoding methods. Ref. Q mentioned that there 
is a maximum 6 dB BER performance penalty for QPSK 
modulation when A = and 7^ 0. To the best of our 
knowledge, no quantitative results and concrete explanation 
have been given for the QPSK case for general A and (p. 

Ref. ifTOl used an amplify-and-forward method in which the 
overlapped signals received at the relay are simply amplified 
and broadcasted back to the end nodes. This method does not 
require synchronization between the end nodes, but introduces 
some inefficiency, because the relay's noise is forwarded along 
with the signal to the end nodes. Ref. 17] investigated systems 
in which symbols are aligned but phases are not. It uses QPSK 
for uplinks, but a higher order constellation map {e.g., 5QAM) 
for downlinks when the uplink phase offset is such that the use 
of QPSK will lead to poor performance. Ref. ifTTl investigated 
OFDM PNC. Because of the long duration of OFDM symbols 
in the frequency domain, the symbols from two end nodes can 
be considered as "almost aligned". However, the phase offset 
is not addressed. 

None of the prior methods attempt to look at the symbol 
and phase offset jointly. Furthermore, they try to "deal with" 
asynchrony, rather than to "exploit" asynchrony to gain per- 
formance advantage. As will be shown in this paper, we can 
exploit symbol asynchrony to gain robustness of the system 
against phase asynchrony. 

IV. A Belief Propagation PNC Decoding Scheme 
(BP-PNC) 

This section presents a PNC decoding scheme based on be- 
lief propagation that deals with symbol and phase asynchronies 
jointly. We refer to the decoding method as BP-PNC. We make 
use of the oversampled symbols in (|4]i to construct a Tanner 
graph ||9l as shown in Fig. |2] In the Tanner graph, Yr denotes 
the evidence nodes, and there are 2N + 1 such nodes; 4* 
denotes the constraint nodes (also known as the compatibility 
or check nodes); and X denotes the source (or variable) nodes. 
We aim to decode the combination (a;^[n], a;B[n]) in X. For 
simplicity, we use x*'^ to denote ^^[j]). Note that 

the Tanner graph has a tree structure. This means that BP 
can find the exact solution for the a posteriori probability 
P{x'^'^\Yr) at each source node. Furthermore, the solution 
can be found after only one iteration of the message-passing 
algorithm (see ||9l). From the decoded P{x'^''^\Yr), we can 
then find the maximum a posteriori probability (MAP) XOR 
value, XR[n] — argmax ^ Pr Yr). 

XR[n\^n,n. X A[n]<Sx B[n]=x R[n] 

In summary, BP can converge quickly and is MAP-optimal 
as far as the BER of XA[n] © XRln] is concerned. Note that 
the MAP optimal scheme is also ML optimal because each of 
the possible symbol values for XA[n] and XB[n] are equally 
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where 



Fig. 2. Tanner graph of the BP-UPNC. 

likely {i.e., Pr {xA[n] = 77^) = Pr {xB[n] = '7-b) = j^. 
VAiVb G X, and x is the alphabet of modulated symbols.)- 

1) BP-PNC Design: Let us now consider QPSK modulation, 
in particular, we define x = {1 + ~1 + jj ~1 ~ Ji 1 ~ i} 
as the symbol set. With reference to (|4}, we have XA[n^ ~ 
a/y/2 and XB[n] = hj^f^, where a, € x- Define p^' = 
P(a;A[[fc/2]] =a/^/2,XB[[A:/2J] =6/V2|yfl[fc]). Note that 
here, p^''' is computed based on 2/i?[fc] only, and not on the 
whole Yr. Also, p^''' is fixed and does not change throughout 
the message passing algorithm in the Tanner graph. P2n-i 
1, 2, N are given as follows: 



P^u for n 



P2n-i = P ( 2^aW = ^,XB[n-l]^ ^ 



27ra2/A 



exp ' 



^/2'^^ ^ ^/2 
yg'^pTi - 1] - Re (a + foe^"^) /V2)^ 
2^VA 



exp • 



(2/]^"[2?i - 1] - Im (a + 6e^*) /v^)' 



2a2/A 



(8) 



(4^[2n]-Re(a 



foe^'^)/V2)' 




that except for the 
P2n-i and has 
and 4 



Note 
of 

possibilities for a 
and last symbols 



exp 



have 

1] = 

. ^ ; (yg°[l]-Re(a/V2)) 

277ctVA 1 25^73 

P{xA[N+l] = Q,XB[n] 

, (yg°[2iV+l]-Rc(be^V72)) 

./^"[2jV+l]-Im(fce^'»/V2)) 
2crVA 



= P{xa[N 



andp2}v+i 



first and last symbols, each 
16 possible combinations (4 
possibilities for h). The first 
4 possibiUties, as follows: 

a/A-T^sW = 0|yH[l]) = 

' (j;^"[l]-Im(a/V2)) " 

25^73 



6/V2|yH[27V+l]) = 



exp 



2j Message Update Rules: Given the evidence node values 
computed by ^ and (|9]l, we now derive the message update 
rules for BP-PNC. Since the Tanner graph (in Fig. |2) for BP- 
PNC has a tree structure, the decoding of the joint probability 
can be done by passing the messages only once on each 
direction of an edge |[T2l . As described below, we could 
consider the right-bound messages first followed by the left- 
bound messages. 

We represent messages on the edges with respect to the 
compatibility nodes 5* in Fig. |2]by Qk and Rk- Specifically, 
Qk (Rk) is the message connected to the right (left) of the A:-th 



R ' 



Oy^ 6y, 0>-; 




(a) 



(b) 



(a) 



(b) 



Fig. 3. Messages passing from left 
to right. 



Fig. 4. Messages passing from 
right to left. 



compatibility node. Pk = (pl^^'^^^ ,pl'^^'~^~^\ ...,pl~^'^~^) 
is a 16 X 1 probability vector, where each component p^''' is 
the joint conditional probability of (a;^[[fc/2]],a;B[[fc/2j]) 
given yR[k] of a particular (a, 6) combination. 
Similarly, - (<il+'-'+^ ,ql^^^-'+\ ...,qlr'^'-^) and 

Rk = {rl+'''+\rl+'--^+\ ...,rl-'-^-') are also 16 x 1 
probability vectors, where g^'' and r^,''' are probabilities 
P {xA[\k/2-]] = a/V2,XB[lk/2\] = b/ V2\yR[l], ...,yR[k - 1]) 
andP{xA[\k/2]]=a/V2,XB[lk/2\]=b/V2\yR[l],...,yR[k]), 
respectively, for the right-bound messages as shown in the Step 
1 (note here, for the left-bound message g^ '' and rj!''' are prob- 
abiUties p(a;rfc/2%Lfe/2j ^ {a,b)/V2\yR[k], ...,yR[2N + 1]) 
and p(xr'=/2l.Lfc/2j ={a,b)/V2\yR[k + l],...,yR[2N+l]) 
correspondingly). 

We omit the index k to avoid cluttering in the following 
discussion of message-update rules. We follow the principles 
and assumptions of the BP algorithm to derive the update 
equations, i.e., the output of a node should be consistent with 
the inputs while adopting a 'sum of product' format of the 
possible input combinations lfT2l . 

Step 1. Update of right-bound messages 

With reference to Fig. |3la), we derive the 
update equations for the right-bound message 

R-^ = (ri+J-i+i,ri+J'-i+^...,ri-J'i-J) based on 
the right-bound message and the fixed message 



-i+j" 



node P = 



,P 



) from the evidence 
). Based on the 



sum-product principle of the BP algorithm, we have 



(10) 



For the input message into the leftmost compatibility node, 
however, we should modify ( fTOl i with 



.a.b 



With reference 
equations 



update 
message R~^ - 



to Fig. 
for the 



P 

we 



, ...,q 



Sb) 
message ( 
based on 



derive the 



., r 



the input 
) of the 



compatibility node. Note that for x = x" " and x' = ", 
the common symbol overlapping in two adjacent samples is 
xsfnl, thus we have 



(11) 



where 6 S x and /3 is a normalization factor to ensure that 



1. If x 



and 



the sum of probabilities ^ g'"''' 

a, hex 

x' — 2;"+^'"+^, then the common symbol is a;^[n + 1]. We 
have 

bex 

where a e x /3 is a normalization factor to ensure that the 
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Fig. 5. Simulations results of BP-PNC with BPSK modulation. 

sum of probabilities ^ q'"^'^ = 1. 

We then update the next evidence node message i?'^ based 
on the newly updated Q'^ in Fig. |3lb) and j/^ in Fig. Oa), 
and so on and so forth, until we reach rightmost node. 

Step 2. Update of left-bound messages 

With reference to Fig. HJa) and (b), we use a similar 
procedure as in Step 1 to update the left-bound messages. 

After steps 1 and 2 above, the values of the messages 
converge, thanks to the tree structure of the Tanner graph. We 
then compute the 4-tuple. 

(\ ^ a, 6 a.h a.b \ ^ a, 6 a.b a.b 

/ , P2n i2n ^2n i 2-^ P2n l2n ^2n ' 

(a,b):a!Sb=l+j {a,b):aeb=-l+j 

Ea,b a.b a.b X ^ a.b a,b a,b | 

P2n i2n ^2n i 2-^ P2n l2n ''2n I ' 

{a,b):a(Sb=-l-j (a,6):a06=l-j / 

where a, 6 e The maximum-likelihood XA[n](BxB[n] is 
given by the a © & that yields the largest element among the 
four elements. Note that we only use P {xAin], XBinWYj^), and 
not P {xA[n + l],XB[n]\YB,), to obtain P {xA[n] XB[n]\YB,), 
because the information required to get P {xa [n] ® xb [n] \Yfj) 
is fully captured in P {xA[n], XB[n]\Yii). 

V. Numerical Results 

This section presents simulation results for BP-PNC. We 
compare the performance of asynchronous PNC with that of 
the perfectly synchronized case [Til, focusing on BPSK and 
QPSK modulations. For comparison purposes, the transmit 
amplitude per symbol in the QPSK case is made to be \/2 
of that in the BPSK case so that they have the same energy 
per bit. In particular, Ei, is the energy per bit in the Fig. |6] 
The ideal synchronous QPSK curve therefore overlaps with the 
synchronous BPSK curve (because when symbols and phase 
are synchronized in QPSK, the system can be considered as 
two parallel BPSK systems). 

A. Summary of results 

Our simulations yield the following findings: 

• For BPSK, the 3 dB BER performance penalty due to 
phase or symbol asynchrony using the decoding methods 
in IHj El is reduced to less than 0.5 dB with our method 
here. 

• For QPSK, the BER performance penalty due to phase 
asynchrony can be as high as 6-7 dB when the symbols 



Fig. 6. Simulations results of BP-PNC with QPSK modulation. 

are aligned, as can be seen from Fig. |6la). However, with 
half symbol misalignment, our BP decoding algorithm can 
reduce the penalty to less than 1 dB, thanks to a "diversity 
and certainty propagation" effect (to be explained). 
• Symbol misalignment makes the system more robust 
against phase asynchrony under our ML decoding method. 

B. Detailed description 

Fig. |5] and Fig. |6] show the simulation results for BPSK 
and QPSK modulations, respectively. The x-axis is the average 
per-bit SNR of both end nodes, and the y-axis is the average 
BER = {BERa + BERb )/2. 

For each data point, we simulate 10,000 packets of 2,048 
bits. We use the synchronous PNC as a benchmark to evaluate 
BP-PNC. As can be seen from Fig. |5ja) and (b), the BER 
performance penalties due to phase and symbol asynchronies 
are less than 0.5 dB for all SNR regimes for BPSK. That is, 
for BPSK, BP-PNC reduces the BER performance penalty from 
3dB in H] to only 0.5 dB. Note that BP-PNC uses complex 
sampling that samples the 2-D information in the In-phase 
and Quadrature-phase components of the signal, whereas the 
scheme in fl] only samples in one dimension. For optimality, 
complex sampling is required in BPSK PNC because of the 
phase offset. 

For QPSK with symbol synchrony but phase asynchrony, 
as can be seen from Fig. HJa), the BER performance penalty 
can be as large as 6 to 7 dB. However, with half symbol 
misalignment, as can be seen from Fig. |6lb), BP-PNC reduces 
the penalty to within 1 dB (compared with the benchmark case 
where symbol and phase are perfectly synchronized). In other 
words, symbol asynchrony can ameliorate the penalty due to 
phase asynchrony. This can be explained by the "diversity and 
certainty propagation" effects elaborated in the next subsection. 
In addition, we note from Fig. |6tb) that when there is a 
half symbol offset, the phase offset effect becomes much less 
significant. Specifically, the spread of SNRs for a fixed BER 
under different phase offsets is less than 0.5 dB. 

C. Diversity and Certainty Propagation 

In QPSK, each symbol has four possible values. Thus, the 
joint symbol from both sources has 16 possible values. The 
constellation map of the joint symbol varies according to the 
phase offset. Fig. [8] shows the constellation map of a joint 
symbol when the phase offset is 7r/4, where the 16 diamonds 
corresponding to the 16 possibilities. For example, a point 
with value 1 + (1 — \/2)j corresponds to the joint symbol 
x"'" = (1 + j, —1 — j) in Fig. |7]due to the phase shift {i.e., 
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Fig. 7. Constellation map for phase asynchronous QPSK PNC (with symbol 
synchi'onization). There ai'e four XORed PNC symbols for QPSK. Constellation 
points of the same color in the figure should be mapped to a same PNC symbol. 
And the amplitude of each symbol is \/2. 

1 + (1 - \/2)j = {1 + j) + (-1 - j)e^''/^)- In PNC, the 
16 possibilities need to be mapped to four XOR possibilities 
for ttie PNC symbol. In Fig. |7] the diamonds are grouped into 
groups of four different colors. The diamonds of the same color 
are to be mapped to the same XOR PNC symbol according to 
XA S) xb- In this mapping process, some of the constellation 
points are more prone to eiTors than other constellation points, 
and the BER is dominated by these bad constellation points. 

With reference to Fig. [T] the eight diamonds within the 
green circle are "bad constellation points". Adjacent points 
among the 8 points are mapped to different XOR values, but 
the distances among two adjacent points are short. By contrast, 
the eight points outside the green circle are "good constellation 
points" because the distances among adjacent points are large. 
When symbols are synchronized {i.e., A = 0), there are 
altogether N joint symbols. On average, half of them will be 
bad constellation points with high BER. 

Now, consider what if there is symbol asynchrony, say A = 
0.5. We have 2N + 1 joint symbols, out of which about N 
will be good constellation points. A symbol from a source is 
combined with two symbols from the other sources to form 
two joint symbols received at the relay. Both the joint symbols 
have to be bad for poor performance. Thus, the diversity itself 
may give some improvement. In addition, there is a certainty 
propagation effect, as explained below. 

Consider a good constellation point, say a;"" = 
(x^M, For this point, we may be able to decode not 

just the XOR, but the individual values of XA[n] and XB[n] 
with high certainty. Now, suppose that the next joint symbol 
^n+i,n _ (^XA[n+l],XB[n]) is a bad constellation point. But 
since we have good certainty about XB[n] from the previous 
good constellation point, the uncertainly in a;"+^'" can be 
reduced by the certainty propagated from x"'" with the BP 
algorithm. In other words, once x b [n] is known, x""^^ " is not a 
bad constellation point anymore. Certainty can propagate along 
successive symbols from left to right, as well as from right to 
left, as shown in Fig. [8] Certainty propagation significantly 
reduces BER. 

VI. Conclusion and Future Work 

This paper has proposed and investigated an optimal max- 
imum likelihood (ML) decoding algorithm based on belief 
propagation, called BP-PNC, for asynchronous PNC. This 
method significantly improves the BER performance of the 
prior methods when there are symbol and phase asynchronies 
in the PNC system. In addition, we find that with symbol 
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Fig. 8. An illustration of certainty propagation. The yellow symbols denote 
the good constellation points and the white symbols with grids denote the bad 
constellation points of Fig. [T] respectively. 

misalignment, BP-PNC desensitizes the BER performance to 
canier-phase offset. In practice, it is relatively easier to control 
symbol timing than carrier-phase offset because of the relative 
time scales involved. Our results suggest that if we could 
control the symbol offset in PNC (e.g., using the method in 
15]), it would actually be advantageous to deliberately introduce 
symbol misalignment so that the system is robust against 
uncontrollable phase offsets. 

Prior to this work, it has often been thought that strict syn- 
chronization is needed for PNC. Our work, however, suggests 
that the penalty due to asynchrony can be nullified to a large 
extent with the BP-PNC decoding method. Our work is based 
on BPSK and QPSK modulations. We believe the same con- 
clusion can be drawn for other modulations, particularly those 
with dense constellations. In addition, channel coding (with 
or without symbol offset in the system) also has the diversity 
and certainly propagation effects of symbol offset expounded in 
Section lV^ Thus, we believe that the incorporation of channel 
coding will also lessen the negative effect of asynchrony in 
PNC systems. These conjectures remain to be investigated in 
future work. 
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